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Abstract 

We propose a new method for reconstructing an effective spin direction of a semi-leptonically 
decayed top quark. The method is simple: for instance, it does not require the spin in¬ 
formation of the antitop quark in a tt event. The reconstructed effective spin is expected 
to be useful in hadron collider experiments. We demonstrate its usefulness in an analysis 
of the top decay vertex. 


After ten years since the discovery of the top quark P], as yet only limited experimental 
data on its properties are available, and details of the nature of the top quark are still to be 
unveiled. So far, there exists no evidence for signihcant deviations from the Standard Model 
(SM) predictions concerning top quark properties. On the other hand, since the mass of the top 
quark approximates the electroweak symmetry-breaking scale, there are hopes that symmetry¬ 
breaking physics may manifest itself through non-standard interactions of the top quark. For 
detailed examinations of top quark interactions, it is expected that the top quark spin can be 
used as a powerful analysis tool. This is because, in the SM (and in many of its extensions), 
the top quark decays before it hadronizes, and the spin information of the top quark is directly 
reflected to the distributions of its decay products. Hence, we may utilize the top quark spin 
for disentangling different top quark interactions efficiently. This is in contrast to the other 
lighter quarks, for which hadronization effects dilute the spin information at the quark level 
severely. 

In a future linear collider experiment, we will be able to utilize the spin of the top 
quark quite efficiently. This is because produced top quarks will naturally be polarized, due to 
parity-violating nature of the interactions in the top quark production process. The top quark 
polarization can even be controlled or raised to high values by tuning the polarization of the 
initial electron beam. 

By contrast, spin reconstruction of top quarks in hadron collider experiments is a non¬ 
trivial task. At Tevatron and LHC, top quarks are produced mainly through tt production 
processes, and these top quarks are known to be hardly polarized j2j. Two types of methods 
have been proposed and studied for utilizing the top quark spin at these colliders. One is to 
take advantage of the correlation between the top quark spin and antitop quark spin in the 
ti events. The other method is to use polarized top quarks produced through the single top 
production process. Unfortunately, so far, not much information has been obtained by applying 
these methods in analyzing real top quark data in the Tevatron experiments, due to intrinsic 
disadvantages of the methods. The only analysis that has been performed is a spin-correlation 
measurement by DO, which put a fairly loose bound on a correlation coefficient jS]. 

Disadvantages in using the top-antitop spin correlation in the ti events are as follows. If we 
analyze the ti events which decay in the dilepton channel (both t and i decay semi-leptonically), 
data statistics is low due to the small branching fraction. Furthermore, there are two missing 
neutrino momenta in each event, which make reconstruction of the event topology non-trivial. 
Instead, if we analyze the ti events which decay in the lepton + jets channel (one of the top 
quarks decays semi-leptonically and the other decays hadronically), we have more statistics, 
but reconstruction of the spin of the hadronically-decayed top quark needs to go through 
complicated procedures. The reconstruction process is affected signihcantly by kinematical 
cuts, and often important information is lost by the cuts. Accurate estimation of the effects of 
kinematical cuts and event reconstruction is crucial as well. The complex procedures bring in 
sizable systematic uncertainties in the top spin reconstruction. 

Disadvantages in using the single top production process are that there are huge background 
cross sections for Wbb and Wbb + jets processes, and that these background cross sections have 
not been estimated accurately. In fact, up to now the single top production process has not 
been observed at Tevatron, as opposed to original expectations, due to lack of data statistics 
and difficulty in the background estimation. 

We would expect that, when a huge top quark sample is available at LHC, these difficulties 
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will eventually be overcome, along with reduction of statistical errors as well as better under¬ 
standings of systematic uncertainties; see e.g. |3]. On the other hand, it is certainly desirable to 
develop another method for top spin reconstruction that can be applied to a high statistics sam¬ 
ple, and that involves small and controlled systematic uncertainties. In this paper, we propose 
a new method that can meet such criteria. This method can be applied to a semi-leptonically 
decayed top quark, without requiring reconstruction of the spin of the antitop quark. We will 
demonstrate usefulness of our method in an analysis of the top-quark decay vertex. A more 
detailed application of our method is given in [H], where sensitivities to anomalous couplings 
in the top decay vertex are studied, using our method, and taking into account realistic ex¬ 
perimental conditions expected at Tevatron and LHC. There, it is shown that our method for 
effective top spin reconstruction is indeed practically useful. 

One may be perplexed, since there is no spin vector (polarization vector) associated with 
an unpolarized state; one may well argue that it is impossible to reconstruct the spin of an 
unpolarized top quark. While this argument is correct on its own, we can still reconstruct an 
effective “spin direction” of a top quark, which is practically useful in analyses of top decays. 

The method we propose is simple and naive. It is based on the following two well-known 
facts: 

(i) Within the SM, the charged lepton in the semi-leptonic decay of a 100% polarized top 
quark is emitted preferentially in the top quark spin direction. In fact, at tree level of the 
SM, the normalized angular distribution of the lepton is given by [HI 

1 dr{t ^ blu) 1-I-COS0 , , 

J\f d cos 0 2 

where 0 denotes the angle between the top spin direction and the lepton direction in 
the top rest frame. A/” represents the normalization constant such that the integral upon 
cos 0 becomes unity. It is known that the one-loop QCD correction hardly modifies the 
above angular distribution [7]. 

(ii) If we include anomalous couplings in the top decay vertex, their effects on the lepton 
angular distribution enter only from quadratic dependences 0 , (All terms linear in the 
anomalous couplings vanish.) Namely, when the anomalous couplings are small, their 
effects are very suppressed. 

Unpolarized top quarks can be interpreted as an admixture, where one-half of them have their 
spins in +n direction and the other half have their spins in —n direction, for an arbitrary 
chosen unit vector n. The directions of the charged leptons from the top quarks with ±n 
spin are emitted preferentially in the ±n direction in the top rest frame. Hence, it seems 
reasonable to project the direction of the lepton pi/\pi\ onto the n-axis and define an effective 
spin direction as sign(fi-p)) x n for each event. According to eq. (Q), in this way we choose 
the correct direction with probability 75% on average. That we can choose any axis n, and 
that any choice is equivalent (if we ignore experimental environment), guarantee the rotational 
invariance of the unpolarized state of the top quark. Due to the above property (ii), the dehned 
direction is hardly affected by the anomalous couplings in the top decay vertex if they are small, 
so that it is appropriate for an effective spin direction. 

Importance of this definition consists in our finding that certain angular distributions of 
the top decay products with respect to the effective spin direction reproduce fairly well the 
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corresponding angular distributions from a truly polarized top quark. This is the case even 
including anomalous couplings in the top decay vertex. This is the main aspect to be addressed 
in this article. 

Provided that produced top quarks are perfectly unpolarized, and provided that we disregard 
effects by kinematical cuts and acceptance corrections, there is no difference on which spin axis 
n we choose to project the direction of the charged lepton. In most part of the paper, we 
consider this ideal case. We will briefly discuss effects of incorporating realistic experimental 
conditions at the end. 

We start by explaining our setup of the top decay vertex including form factors. We assume 
that deviations of the top decay form factors from the tree-level SM values are small. Then 
we consider only those form factors which induce deviations of the differential distributions of 
top decay products at the hrst order in the anomalous form factors. That leaves only two form 
factors and in the limit mi, —> 0 and for onshell W, although the most general tbW 
vertex includes six independent form factors UD]: 


Kr,t = -^V„,U(P,) 


ft pl-' 4?^ f? 
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( 2 ) 


where = (1 =f 75)72 are the left-handed/right-handed projection operators; k is the mo¬ 
mentum of W. For simplicity, we further assume that and are real.* At tree level of the 
SM, f[ = 1 and = 0. We will be concerned only with the top decay process t bW, where 
the value is fixed, therefore, we treat the form factors as constants (couplings) henceforth. 

Using the above decay vertex and taking the narrow width limit for W, the differential 
decay distribution of W and I in the semi-leptonic decay of a 100% polarized top quark is given 
by [10] 
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Here, is the Fermi constant. 9w is defined as the angle between the top polarization vector 
and the direction of W in the top quark rest frame. 9i is defined as the lepton helicity angle, 
which is the angle of the charged lepton in the rest frame of W with respect to the original 
direction of the travel of W. (pi is defined as the azimuthal angle of I around the original direction 
of the travel of W. A schematic view of the angle definitions is shown in Fig. ^ The above 
differential distribution contains fully differential information on the decay t bW —> blu. 

The corresponding differential distribution with respect to the effective spin direction can 
be computed in the following way. An arbitrary unit vector n is chosen as the spin axis in the 


*We note that the anomalous couplings for the right-handed bottom quark and the CP-odd anomalous 
couplings, which we neglect here, are severely constrained indirectly from the measurements of Zbb vertex at 
LEP and of 6 —> 57 process [O]. This may provide another justification for neglecting these form factors in our 
simplified analysis. 
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Figure 1: Schematic view of the angle definitions. 


top rest frame. We denote the charged lepton momentum as pi in the same frame. Then, if 
n ■ pi > 0, we dehne the effective spin vector to be Seff = n, whereas if n ■ p/ < 0, we dehne the 
effective spin vector to be Seg = —n. The angle ©eg between Sgg and pi is given by 


cos ©eff = 


^eff ' Pi 
\pi\ 


1 + Pw cos 9i 


^sin 9i cos pi sin 9w + 


cos 9i + Pw 

\/l ~ Pw 



(5) 


where Pw = {'rnp — M^)/{7rp + M^) denotes the velocity of W in the top rest frame. 9w and 
pi are dehned as in Fig. [T]with respect to Seg (instead of the top polarization vector)."^ Thus, 

'dT{t blv) 

dcos 9wd cos 9idpi 


pff. snin 


dTjt -^bW^ blv) 
d cos 9wd cos 9idpi 


X 2 6* (cos ©eg). 


( 6 ) 


_ iinnol. 


9{x) represents the unit step function. Here, the decay distribution from an unpolarized top 
quark is given by 


dT{t ^bW^ blv) 
d cos 9wd cos 9idpi 


unpol. 
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Obviously, it is independent of 9w and pi, since there is no reference spin vector. Therefore, 
the dependences on 9w and pi of the differential distribution with respect to the effective spin 
direction enter only through the step function 6*(cos©eg) on the right-hand-side of eq. (jHl). 

At this fully differential level, [dT/dcos9wd cos9idpi\ef[ spi^ [eq. Q] is only a crude approx¬ 
imation to dV/dcos9wdcos9idpi [eq. Q]. It can be seen, for instance, from the existence of 
the step function or from the factorized form of the dependences on {fp, f^-) and on {9w, Pi) 
in eq. neither of which is in the structure of eq. dSD. 

Let us integrate over pi and compare the double angular distributions with respect to the 
true and effective spin directions: 

Ht would be more accurate to denote these angles as 0w,eS and pi^eS, but to avoid illegibility we use the 
same notation as in eq. IpHl. 
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Figure 2: Normalized double angular distributions for (/i, / 2 ) = (1,0) (a) using the true spin direction 
and (b) using the effective spin direction, corresponding to eqs. (|S1) and Q, respectively. They are 
normalized to unity upon integration. 
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COS 9i + f3w cot 9w 
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0 if a; > 1 

27r if a; < —1 . (10) 

71 — 2 arcsin x if—l<x<l 


Numerically these two distributions become reasonably close to each other.1 This is demon¬ 
strated in Figs.[21(a)(b), in which both double angular distributions are displayed for (//', f^) = 

(1,0) (tree-level SM). The distributions are normalized to unity upon integration. Qualitative 
features of the bulk distribution shape of dT/d cos 9wd cos 9i are reproduced by [dF /d cos 9wd cos dijeg. spin- 


It has been known that the double angular distribution dT/d cos 9wd cos 9i is useful for 
probing the anomalous coupling HH. To see sensitivities to the anomalous couplings semi- 

Tt is not obvious from the explicit formulas for the distributions. In particular, eq. 01 still has a factorized 
form concerning the dependences on (/f, f^) and 9 w, which is different from eq. 0. 
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quantitatively, we divide the phase space into four regions as 


Region A : 

-1 

< 

cos 9w 

< 

0 

and 

-1 

< 

cos 9i 

< 

0 

Region B : 

-1 

< 

cos 9w 

< 

0 

and 

0 

< 

cos 9i 

< 

1 

Region C : 

0 

< 

cos 9w 

< 

1 

and 

-1 

< 

cos 9i 

< 

0 

Region D : 

0 

< 

cos 9w 

< 

1 

and 

0 

< 

cos 6 *; 

< 

1 


and define the event fraction in each region by 

Ri=M-^ ' 

where 
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(13) 


represents the top-quark partial width to blu. Each Ri is a function of f^/fi- We also define 
the event fractions Rf^ in the same manner using the effective spin direction instead of the true 
spin direction. 

We compare the dependences of Ri and Rf^ on f^/fi in Figs. |Sl^a)(b). From the figures, 
we see that major features of the f 2 /fi dependences of Ri are reproduced by In fact, 
the dependences of are consistent with the observation that, if we use the effective spin 
direction, we misidentify the correct spin direction with 25% probability on average. Namely, we 
misidentify Region A with C, and Region B with D, so if we combine RiS in Fig. Ha) reweighting 
them with this misidentihcation probability, we obtain the curves similar to those plotted in 
Fig. Hb). Since the f^lfi dependence of the most sensitive event fraction is about half 
of that of Ra, if we use the effective spin direction, we would expect a sensitivity to f^/fi 
roughly half of what would be obtained with the true spin direction. A closer examination of 
the sensitivities to f^/fi incorporating realistic experimental conditions is given in |H]. 

The one-loop QCD correction to dV/d cos 9wd cos 9i has been computed in na, A large 
part of the correction goes to a variation of the normalization of the partial decay width, 
which amounts to about 9%. On the other hand, the correction to the normalized double 
angular distribution is at the level of 1-2% or less. Although the one-loop QCD correction 
to [dr/dcos^vi/dcos^/jeff, spin has not been computed yet, we expect that it would not be very 
different from the correction to dV/d cos 9y/d cos 9i. If this is so, we may be able to measure the 
QCD correction to the normalized double angular distribution at LHC, provided that a good 
understanding of systematic uncertainties is possible; cf. j5]. 

We may also compare the angular distributions dV/dcos9i and [dT/dcos^jjeff.spm, where 
9i denotes the angle between the direction of particle i and the top polarization vector or the 
effective spin direction SeS in the top rest frame. The normalized angular distributions are 
shown in Fig. 0] for i = b,W,i> and = (1,0). (The lepton angular distributions are 

trivial, so we do not show them here.) The angular distributions with respect to the true spin 
direction depend linearly on cos 9i. In this case, it is customary to parametrize a normalized 
angular distribution by |(1 -|- a* cos 6 *j) and refer to as a correlation coefficient. Since b and 
W are emitted back-to-back in the top rest frame, ab = —aw The correlation coefficients 
ttj corresponding to the true spin direction have been computed in HB]. On the other hand. 
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Figure 3: Event fractions in the four regions as functions of (^) using the true spin direction 

[Ri) and (b) using the effective spin direction (iZf®). 


the angular distributions with respect to Seff are not given by linear functions of cos6i. Their 
analytic expressions are complicated, which we do not present here. Numerically, the angular 
distributions for b and W are close to linear shape, while that of u is considerably different 
from linear shape close to cosOi^ = ±1. 

If we approximate [dF /d cos 6*j]eff. spin for i = b,W hj linear functions of cos 9i, the correlation 
coefficients \ab\ and \aw\ for the effective spin direction are about twice larger than those for 
the true top spin direction. As for the angular distribution of v with respect to Seg, on average 
the slope of the distribution is steeper (the correlation between the neutrino direction and Seff 
is stronger) than that of dV/dcosdy. These enhancements in the angular correlations, if we 
use the effective spin direction instead of the true spin direction, stem from purely kinematical 
origins. It can be understood as follows. Consider a hypothetical case, in which no correlation 
between the true spin direction and direction of W exists (the decay is isotropic). Even in 
this case, there is a positive correlation between the effective spin direction and W in the top 
rest frame, since the charged lepton is emitted more in the direction of W due to the boost 
by W. Similarly, (hypothetically) even in the absence of any correlation between the true spin 
direction and neutrino direction, there is a negative correlation between the lepton direction 
and neutrino direction in the top rest frame, since they are 100% anticorrelated (back-to-back) 
in the W rest frame. These kinematical effects bias the angular correlations to be stronger if 
we use the effective spin direction. 

We also examined the /^//f dependences of the angular distributions dV/dcos6i and 
[dF/d cos spin- The f^/fi dependences of the latter distributions are much weaker than 

those of the former distributions. The dependences of [dF/d cos djjefr. spin for i = b,W are 

reduced as compared to the f^/fi dependences of the double angular distribution eq. 0- This 
is due to cancellations of /^//Z dependences between and and between and 
see Fig. infb). Insensitivity of [dF/d cos 6*,^]eg. spin to f^/fi stems from a strong (anti)correlation 
between the effective spin direction (or the lepton direction) and the v direction. Since the 
f^/fi dependences of [dF/d cos djjeg. spin are weak, it is much more advantageous to use the 
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Figure 4: Normalized angular distributions for b,W,u in the top rest frame for (//", f^) = (1, 0), using 
the true spin direction (solid lines) and using the effective spin directions (dashed lines). 


double angular distributions [eq. (jHI)] or event fractions Rf for gaining sensitivities to 

Up to now, in defining the effective spin direction, we assumed that the initial top quark 
is completely unpolarized and neglected experimental environment. In practice, under realistic 
experimental conditions, different choices of spin basis (axis) n lead to different distributions 
of decay products. Effects of kinematical cuts are by far the largest. Based on detailed Monte 
Carlo simulation studies incorporating realistic experimental conditions expected at Tevatron 
and LHC, it is found that the top helicity axis pt/\pt\ defined in the ti c.m. frame (opposite of 
the direction of i in the top rest frame) is an optimal choice for the spin axis n. Other choices, 
such as beamline axis and the off-diagonal spin basis jEj, turn out to be inappropriate, since 
the original distributions [e.g. Fig. |2l(b)] are strongly distorted by the effects of kinematical 
cuts, and also because the sensitivities to the anomalous couplings are much reduced. This 
can be understood as follows. If we choose the beamline axis, the small Et (transverse energy) 
and large |? 7 | (pseudorapidity) regions correspond to the regions cos9w — ±1, and events 
that fall into these kinematical regions are rejected by cuts such as the requirements for the 
minimum transeverse energy {E^ cut) or acceptance correction (Ir^l cut) for the lepton and jets. 
In particular, events in the kinematical regions most sensitive to a variation of close 

to {cos 9w, cos 9i) = (1,1) and (-1,-1), are lost. At Tevatron, the status of the off-diagonal 
basis is similar to the beamline axis, since the off-diagonal basis is not very different from the 
beamline axis. (At LHC, there is no good definition of the off-diagonal basis.) On the other 
hand, if we choose the top helicity axis, after integrating over all top quark directions, effects of 
the cuts are averaged over and no signihcant distortion from the original distribution is found. 
See jHj for details. 

As is clear from the above dehnition, our method can be applied not only to the hadron 
collider experiments but also to a future e+e“ collider experiment. Nevertheless, the primary 
motivation of our proposal is to use this method at the current Tevatron experiment and at 














LHC. 

In summary, we proposed to reconstruct an effective spin direction of a semi-leptonically 
decayed top quark as the projection of the lepton direction onto an arbitrary chosen axis in 
the top rest frame. The reconstruction method is simple so that it would be feasible in hadron 
collider experiments. We demonstrated that this spin direction can be used to probe anomalous 
couplings in the top decay vertex, through measurements of a double angular distribution or 
event fractions Under realistic experimental conditions, the top helicity axis seems to be 
an optimal choice for the spin axis. 
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